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Abstract 

. The wavefunctions in phase-space representation can be expressed as entire functions 
of their zeros if the phase space is compact. These zeros seem to hide a lot of relevant and 
explicit information about the underlying classical dynamics. Besides, an understanding of 
their statistical properties may prove useful in the analytical calculations of the wavefunc- 
tions in quantum chaotic systems. This motivates us to persue the present study which by 
a numerical statistical analysis shows that both long as well as the short range correlations 
exist between zeros; while the latter turn out to be universal and parametric-independent, 
the former seem to be system dependent and are significantly affected by various param- 
eters i.e symmetry, localization etc. Furthermore, for the delocalized quantum dynamics, 
the distribution of these zeros seem to mimick that of the zeros of the random functions 
as well as random polynomials. 
PACS numbers: 05.45 +b, 03.65 sq, 05.40 +j. 
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1. Introduction 

In this paper, our aim is to study the statistical distribution of zeros of "quantum 
chaotic" wavefunctions, expressed in phase-space representation. The wavefunctions of 
quantum chaotic systems have so far remained a relatively less explored (as compared to 
spectra) area notwithstanding their importance in quantum mechanical as well as semi- 
classical analysis. The few of their known properties are the microcanonical nature of 
wavefunctions if the underlying classical dynamics is chaotic [1], the scarring behaviour 
[2,3], dynamical localization [4] or a few statistical (numerical) studies of the components 
of these wavefunctions [5,6]. The relevant information hidden in the semiclassical wave- 
function about the underlying quantum-classical correspondence is not yet fully extracted. 
This motivates us to pursue the present problem. As intution suggests us, this can be best 
studied by analyzing various properties of the wavefunctions in the phase-space where 
quantum dynamics in limit h — > can be compared with the underlying classical dynam- 
ics. One such analysis involves the study of the distribution of zeros (or nodal patterns) 
of wavefunctions (in phase space representation). 

The interest in the phase space study of the nodal patterns arises following a study by 
Leboeuf and Voros [7,8] where they show that the coherent state representation of a wave- 
function (i.e Bargmann or Husimi function) has a finite number N of zeros in phase space 
(N being the dimension of Hilbert space) and the state of the system can be completely 
expressed in terms of these zeros. Furthermore, their numerical studies on quantum maps 
(e.g Baker and kicked rotor) have also revealed that the zeros, in the semiclassical regime 
h — > 0, condense on lines for classically integrable systems whereas for strongly chaotic 
quantum maps they appear to diffuse fairly uniformly over the phase space; this behaviour 
can be explained semiclassically [7,8]. The only exception to this behaviour, so far, has 
been seen for one eigenfunction in N = 2 k case of the Baker map (strongly chaotic) where 
zeros have a dominantly linear distribution. In general, the distribution of zeros seems to 
mimic the underlying classical dynamics in semiclassical limit (e.g. for those chaotic cases 
of the kicked rotor, where quantum dynamics is localized in momentum space, zeros were 
observed to be localized). This makes it relevant to seek further information about their 
behaviour. 

In principal, we have an analytical formulation available to write the eigenfunction 
directly in nodal pattern representation, (a polynomial for spherical phase space and a 
complicated function for torus and cylindrical phase space) [7] . But this formulation does 
not help us very much as it depends on the knowledge of the position of each zero in 
phase space which is analytically difficult to evaluate. Still it can be used to determine the 
statistical behaviour of the eigenf unctions if the statistical properties of zeros are known. 

The random nature of chaotic wavefunction and its polynomial form in Bargmann 
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representation (in a spherical phase space only) has motivated people to seek analogy 
between the distribution of it's zeros and roots of random polynomials (later referred as 
RP) [9]; recently the statistical properties of the latter case were analytically explored by 
us [10]. Our results showed an excellent agreement with the corresponidng ones for the 
kicked top and the kicked rotor. It was a little surprising in the latter case where dynamics 
is confined to a torus and an arbitrary wavefuction is a complicated function instead of a 
polynomial; one should expect the zeros of random functions (superposition of functions 
with random weight, later referred as RF) to be a better model. This seemed to hint that 
the zeros, while distributing themselves in semiclassical limit, do not distinguish between 
spherical and torus phase space [10]. In other words, the zeros of RF and RP seem to have 
the same distribution properties. As we show in this paper, this is indeed the case. We 
also extend the above study to understand the effect of various parameters (i.e symmetry- 
breaking, localization etc.) on the distribution of zeros. We expect this due to our prior 
knowledge of the significant influence of these parameters on quantum dynamics [4,11,12]. 

We choose the kicked rotor system for this purpose as it has been an active model 
of research, containing a variety of features such as localization, resonance, dependence 
of the spectra on number theoretical properties etc and has been used as a model for a 
very wide range of physical systems [4]. Besides, the kicked rotor can also display discrete 
symmetries such as parity and time reversal [5,11,12]. This will enable us to study the way 
these symmetries reveal themselves in the distribution of zeros. 

This paper is organized as follows. The section 2 contains a brief review of quantum 
properties of kicked rotor. The section 3 briefly describes the various phase space repre- 
sentations and their advantage over each other. The section 4 deals with the numerical 
analysis of distribution of zeros for kicked rotor for various parametric conditions, followed 
by a conclusion in section 5. 

2. The Kicked Rotor: Classical and Quantum Dynamics 

The kicked rotor can be described as a pendulum subjected to periodic kicks (with 
period T) with the following Hamiltonian 

H= (P+ 2 7) +^- cos(2n6 + 6 ) £ 5(t-nT). (1) 

n= — oo 

where K is the stochasticity parameter. The parameter 7 and 6>o are introduced in the 
hamiltonian in order to mimick the effects of the time-reversal (T) and the parity (P) 
symmetry breaking in the quantum hamiltonian while acting in a finite Hilbert space. 

Floquet's theorem enables us to describe the related quantum dynamics by a discrete 
time evolution operator U = GBG where B = exp(— iK cos(27r6> + 9q)/4ix 2 Ti) and G = 
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exp(— i(p + r y) 2 /4h). Here 9 and p are position and momentum operator respectively, with 
p having discrete eigenvalues p\l >= l%\l > due to periodicity of 9, (9 = 9 + 2ir). The 
nature of the quantum dynamics and therefore the statistical properties of the associated 
quantum operators depend on h and K. The choice of a rational value for hT/2ir results 
in a periodicity also for the momentum space (p = p + 27rM or / = l + N) and, therefore, in 
discrete eigenvalues for 9 (9\n >= ^f^ln >). The quantum dynamics can then be confined 
to a 2-dimensional torus (with finite Hilbert space of dimension N) and U can be reduced 
to a finite iV-dimensional matrix of the form [7,4] 

3 

exp [-i (n 2 h(l 2 + I' 2 ) - tt 7 (/ + /'))] exp 

where N\ = {N — l)/2 (with N odd) and m, n = -N u -N x + 1, N x . 

In contrast to the classical dynamics, the quantum dynamics can be significantly 
affected by the relative values of three parameters, namely, K, % and N. It has already 
been pointed out [4] that, for K 2 > Nh 2 (= 2nM/T), the eigenstates are fully extended 
in momentum space (later referred as delocalized case); it is also termed as strong chaos 
limit as here the kick is strong enough to make the classical motion ergodic in whole of 
the phase space [12]). We also know that the spectrum as well as the distribution of 
eigenvector-components in this case, with parameters 7 and 6>o chosen to preserve either 
exact or partially violated symmetry, can be modeled by the random matrix theories 
(RMT) [5,6,11,12]. Under this limit, the quantum dynamics has a time reversal symmetry 
T for 7 = and a parity symmetry P for #0 = 0. Though the T-symmetry may be 
violated for 7 7^ 0, still a more generalized antiunitary symmetry S = TP = PT can be 
preserved in the system if #0 = [5,11]. By a slow variation of these parameters, one can 
realize the various intermediate stages of the statistical properties of quantum operators, 
for example, the Poisson — * COE transition of the spectrum for K-variation, COE — > CUE 
for 7-variation with 6>o fixed at nonzero value etc. In the opposite limit of weak chaos, 
namely, K 2 « NU 2 (corresponding to a diffusively covered phase space with a non- 
uniform distribution of periodic orbits), the eigenstates localize in the momentum space 
and one obtains a Poisson distribution for the spectrum (later referred as localized case). 
Moreover, here the symmetry breaking may not be realized in the quantum dynamics just 
by making 9q or 7 sufficiently non-zero. 

3. The Phase Space Representations of Wave Function 
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The phase space representation of the eigenfunctions of quantum mechanical operators 
can give excellent information about the semiclassical description of a quantum state. A 
semiclassical wavefunction can be constructed, first by writing the quantum state as a 
function of the parameter h, and then by exploiting the fact that this function should 
approach the corresponding classical function in the limit fi — > [13]. This is easy to be 
observed in the phase space representation of quantum state, where the quantum dynamics 
(the Heisenberg equation), analyzed in terms of density operator p = \ip >< explicitly 
appears as a deformation of classical dynamics (the Liouville equation). This enables one 
to compare the quantum function with the corresponding classical function, order by order 
in h. The inversion W — > if} then gives, in principle, the quantum state for various orders 
in fi. The two most widely used phase space techniques are the Wigner representation W w 
defined as 

W w (q,P, h) = (2nh)~ N J - r -Mq + ^)e^ h d N r, (4) 
and the Husimi representation Wh given as 

W h (q,p,h)= l<Z ^ >l \ (5) 
where \z > is a plane coherent state described by following standard formulas, 

|*>=e*» + |0>, a+ = ii^, z= q -^, (6) 

with q and p as the position and momentum eigenvalues of the operator q and p respectively 
and its overlap with position wave function given as follows, 

(q\z) = (Trh)- N '*exp[(-(z 2 + q 2 ) + V2zq) (7) 

Both of these representations have a nice semiclassical behaviour, converging to the classi- 
cal phase space distribution in the limit h — > 0. But the determination of this convergent 
form of W w and Wh can be done only after smoothing out the oscillations (stronger in 
the case of W w ) existing in both cases. In general, however (i.e for chaotic cases) even 
this knowledge is not sufficient to reconstruct the wavefunction itself by inverting the map 
W — > ip. The reason is that the quantum phase, which controls the fast oscillations, is 
obliterated in the limiting process and can not be easily regenerated. Moreover these rep- 
resentations also run into difficulties when applied to compact phase spaces. It is because 
the existence and uniqueness of the construction of Wigner representation, which imposes 
itself by symmetry on a linear phase space, is no longer clear on a compact spaces [7]. 
One does not face these difficulties with Husimi representation which only needs coherent 
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states. But the usual plane coherent state, used in Husimi formulation, becomes nonana- 
lytic when restricted to a compact space [7]. Thus, for our purpose, we should choose a 
coherent state with analytical properties in compact space. 

As discussed in [7], the analyticity of quantum state can be preserved by expressing 
it in "Bargmann representation". This representation diagonalizes the creation operator 
z and expresses a state as function ip(z) =< z\ip > with z + = h-^ (similar to Schorindger 
representation ip(q) of quantum state, obtained by diagonalizing the position operator q 
with p = —ih-^). The function ip(z), describing the overlap of wavefunction ip with the 
coherent state \z >, belongs to a Hilbert space of entire (analytic) functions with finite 
norm given as follows, 

< i){z)\i){z) >= (27rh)- 1 J \ifj(z)\ 2 e-W 2 / n dqdp. (8) 

Furthermore the Husimi function can be expressed in terms of ip(z) as Wh = \ < z\ip > 
| 2 / < z\z >. Thus one can consider ip(z) as a sort of phase space representation for the 
wave vector >. The advantage of working with tp(z) is that it preserves the information 
about the quantum phase explicitly. 

It has been shown in ref.[7] that the antianalytic coherent states for compact space 
such as the torus (the case considered here) can be obtained by applying periodic boundary 
conditions on the antianalytic coherent states \z > on the plane. Thus the imposition of 
the appropriate periodic boundary conditions on the kernel (5), depending on the nature 
of the compact space, leads to a Hilbert space of entire functions ip(z). For our purpose, 
we only mention here the case of the torus, where application of the periodic boundary 
conditions in both the q as well as p directions and use of the condition 2t\NU = 1 results 
in following form for ij)(z) 

N-l 

^( z ) = Yl ^(*») « z \ qn ( 9 ) 

n=0 

where ipnS are the eigenfunctions of operator U in the coordinate representation and 
<< z\q n » is projection of the coordinate eigenfunction \q n > in the z-space, under the 
periodic boundary conditions. 



« z\q n »= (2iV) 1/4 exp 



-2nN((z 2 +q 2 n )/2-V2zq r 
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(10) 

Here coefficients ip(q n ys are complex under no-symmetry condition. But presence of a T- 
symmetry makes them real as now if ip(z p ) = then ip(Tz p )=0 too [10]. The existence of 
an additional symmetry such as parity imposes further restrictions on ^>(g n )'s; ip(q n ) = 
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if n is odd and ip(q n ) 7^ if n is even. Further, by using the periodic boundary conditions 
leading to this compact space, it can be shown [7] that every function ip(z) has exactly N 
zeros in the fundamental domain (defined by 5q = 1, Sp = 1). 

It is worth reminding ourselves here that, for a spherical phase space, << z\q n » 
takes the form of a polynomial [7] and therefore the choice of ^(qvj's as Gaussian random 
variables leaves i(j(z) as a random polynomial. In ref.[10], we had analytically calculated 
the joint probability density of these zeros as well as the correlation function and compared 
the result with numerical results for kicked top (which has a spherical phase space); we 
found a good agreement. But as obvious from eq.(10), for a torus phase space, ip{z) is not 
a polynomial. A choice of ^(g n )'s as random variable gives rise to a random function (later 
referred as RFC for complex and random coefficients, RFR for real and random coefficients 
and RFA if the coefficients are real, random and non-zero for even n and zero for odd n). 

4. Numerical Analysis 

With this knowledge, we can now proceed in following steps, (i) numerical compu- 
tation of eigenf unctions by diagonalizing the matrix U (eq.(3)), (ii) calculation of their 
Bargmann transforms, using eq.(9), (iii) numerical computation of the position of each 
zero using Cauchy's integral (see [7]), (iv) the, so obtained, A 2 zeros (A for each of the 
A eigenfunctions) are then used for evaluation of the various statistical measures. For the 
reasons explained above, We repeat the steps (iii) and (iv) also for random functions cases, 
with normally distributed coefficients ip(q n ys. 

To see the effect of various parameters on the statistical behaviour of zeros, the above 
steps are repeated for many of their combinations. Here, untill otherwise stated, we always 
take (K = 20000, M = 1, A = 199,7 = 0.7071, O = n/2N) for delocalized case with no- 
symetry, (K = 20000, M = 1, A = 199, 7 = 0.0, 6 = tt/2A) for delocalized case with one 
symmetry (time-reversal), (K = 20000, M = 1, A = 199, 7 = 0.0, 9 = 0.0) for delocalized 
case with two symmetries (namely time-reversal and parity) and (K = 10, M = 30, A = 
199, 7 = 0.7071, #0 = ?r/2A) for localized case with no symmetry. To improve the statistics, 
we superimpose the statistics for N zeros of each of the N eigenfunctions and set T=l. 

We start by computing the density p of zeros, distributed over the phase plane of unit 
area. We know that the eigenfunctions of RM ensembles, in Bargmann representation, have 
the form of the random polynomials. The nature of the coefficients in these polynomials 
depends on the symmetries of the ensembles (e.g complex for no anti-unitary symmetry 
and real for 1 anti- unitary symmetry). Recently it was shown that the zeros for RP 
with complex coefficients are uniformly distributed in the phase space for any arbitrary N 
while, in case with real coefficients, this uniformity no longer survives due to an increased 
concentration of zeros on real axis [10]. The same should hold good for the eigenfunctions 
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of quantum systems with strongly chaotic classical dynamics where a typical wavefunction 
randomly goes all over the phase space and has a uniform intensity distribution except for 
a few scars (high intensity regions) [1,2]. This expectation has its roots in the success of 
RMT in modelling "quantum chaotic" spectra; we hope the same for "quantum chaotic" 
eigenfunctions too. In fact, the figure 1(a) which shows the distribution of 99x99 zeros of all 
99 eigenfunctions of QKR (K = 20000, M = 1, N = 99, 7 = 0.7071, O = n/2N) strongly 
suggest the uniformity of density in the phase plane. The various parameters here are 
chosen so as to ensure the underlying classical dynamics to be strongly chaotic and removal 
of all symmetries from the quantum case [5,11,12]. The figure 1(6) shows the case with 
a chaotic classical dynamics while the quantum dynamics is localized in the momentum 
space (K = 10, M = 30, N = 99,7 = 0.7071, 6 = tv/2N). Here again the zeros seem to 
be randomly and uniformly distributed, except for an increased concentration on parity 
symmetry axis q = 0.5 (reasons explained later). The apparent randomness and uniformity 
of the distribution of zeros (away from the axis) in this case is indeed misleading and occurs 
due to the averaging effects produced by superposition of many wavefunctions. In fact, the 
figure 2(a), where the Husimi distribution of a single eigenfunction is plotted, shows that 
the zeros for the localized case, in general, form patterns; for few of the eigenfunctions, 
zeros seem to be localized too. It seems that the maximum intensity regions, in order to 
concentrate in some part of phase space, compel minimum intensity regions to distribute 
in an ordered way. This pattern formation of zeros is more obvious for the eigenfunctions 
of the localized dynamics with a symmetry (figures 2(6), 2(c)). The figure 2(d) shows 
the Husimi plot of an eigenfunction for delocalized case. Here zeros of even a single 
eigenfunction seem to be uniformly distributed. In fact, the distribution approaches the 
uniformity in N — > 00 i.e the semiclassical limit (unlike the RMT case where uniformity 
exists for all N) and therefore mimicks the behaviour of underlying strongly chaotic classical 
dynamics. 

In order to see whether the density p(p, q) is separable in variables p and q or not, we 
calculate the coefficient of correlation C(m,n), defined as follows, 

n( , < P m q n >-<p m ><q n > m , 
C(m,n) = - — — — r , (11) 

where cr 2 (r m ) is the variance of mth moment of the variable r. To see the effect of 
various parameters as well as dimensionality, we calculate C(m,n), (m + n < 8), for 
delocalized case with no symmetry, one symmetry and two symmetry respectively and 
also for localized case with no symmetry, with each case studied for four values of N, 
namely, N = 61,99, 149, 199. We find that, for all the parametric cases, the correlations 
between p and q variable decrease with N. In fact, the correlations for lower moments are 
much less than one sample error even for N = 61 (with C(m, n) < 10~ 5 for m + n < 5 
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and sample error =10~ 2 ), thus indicating their near-absence even for small iV-values. The 
correlations for higher moments (m + n > 5) are although not negligible for small N- 
values but decrease rapidly with increasing N and can safely be assumed to be zero in 
limit N -> oo, (e.g C(4,4) = 0.018 ± 0.016 for AT = 61 and 0.003 ± 0.005 for N = 199, 
for delocalized case without symmetry). This implies that, for the delocalized as well 
as the localized case, we can study the distribution of zeros in either p or q direction 
without taking care of q or p variable, respectively. Thus, for a check on uniformity, it is 
sufficient to study the measure < p > q and < p > p (where <> r implies the average over 
variable r). As shown in figures 3(a) and 3(6), the averaged density in both q (averaged 
over p) and p directions (averaged over q) for delocalized case is nearly a constant, thus 
indicating a uniform distribution while, for the localized case, it oscillates strongly about 
an average value. Here the increased average density around some values of p and q implies 
the preference of zeros for this value; this behaviour may be attributed to the localization 
and pattern formation of zeros. For the delocalized case, the oscillation strength decreases 
rapidly with increase in N which indicates the uniform density also in p-direction, in large- 
N limit. For comparison, we have also plotted the corresponding RF case in both figures 
3(a) and 3(6); the set of normally distributed ip(q n ys is generated bu using "GASDEV" 
subroutine [14]. (note, the similarity of the delocalized case with the RF case in figure 
3(a),(b)). 

To see the effect of symmetry on the distribution of zeros, we study it for various values 
of parameters 7 and 9q. The figures 4(a) and 4(6) show the distribution of the 61 x 61 
zeros of all 61 eigenfunctions of QKR for (K = 10, M = 30, N = 61, 7 = 0.0, 9 = 0) and 
the case (K = 10, M = 30, N = 61,7 = 0.0, 9 = n/2N). Here the first set of values 
correspond to preservation of both T and P symmetries in the system whereas in the 
second case only T symmetry is preserved. The figure 4(a) shows the increased density 
of zeros on the two symmetry axis (namely, p = 0.5 and q = 0.5) while in figure 4(6) 
the zeros are concentrated only on one symmetry axis; the concentration of zeros also on 
p = 1.0 is due to the torus boundary conditions. Moreover zeros on either side of the 
symmetry axis seem to avoid it, thus creating a hole close to and around the symmetry 
axis. This phenomenon, which occurs in the case of RP (and RF) too, can be explained 
as follows: if ip(z p ) = 0, then ifj(Rz p ) = 0, where R refers to a particular symmetry in the 
quantum system and Rz p the .R-symmetry counterpart of z p . This implies that, given z p 
a zero of wavefunction -0, Rz p will also be its zero; zeros either come by pairs symmetric 
with respect to the symmetry-axis or they are single and lie on the axis. Although this 
destroys the uniformity of distribution of zeros on a global level, locally and away from 
the symmetry axis the uniformity is still preserved (figures 3(c, d), 4(o, 6)). It should also 
be noticed in figure 1(6) that though both 7 and 6>o are non-zero still the density of zeros 
is higher along q=0.5 line while this behaviour is absent in delocalized case for the same 



9 



values of symmetry breaking parameters (figure 1(a)). It seems that for localized case, a 
more generalized symmetry remains preserved even for non-zero values of 7 and 9q. 

In order to analyze the correlations between these zeros as well as to know the extent 
of randomness, we compute the 2-point correlation function R 2 (r) . It is defined as the joint 
probability-density of finding a zero at a radial distance x from the origin and another at 
x + r, averaged over all x and 9 where 6 is the angular dependence of the zero. Here 
r is measured in units of averaged spacing and i?2( r ) is unfolded in such a way so that 
average number of zeros at a distance r from a given zero does not depend on r. In 
Ref.[10], we had analytically obtained R-zir) for zeros of random polynomials; a numerical 
comparison of this formulation with that for kicked spin system as well as kicked rotor 
showed an excellent agreement. Here we numerically compare Rzir) for delocalized-QKR 
with that for the random functions (with complex coefficients to model the no-symmetry 
case and real coefficients for time-reversal symmetry preserving case) and find the latter 
to be a good model for the former; see figures 5(a),(b),(c). The agreement of results 
also confirms that zeros for both random polynomials and random functions, in N — > 00, 
have the similar correlations at least up to the second order. The histograms in figures 5 
(a),(b),(c),(d) also indicate that, for short distances, the two-point correlation increases as 
a power law with distance, attains a maximum value and then oscillates rapidly around 
this average value when distances are large. Thus, for short distances (r < 1), zeros 
seem to repel each other ; the repulsion is weaker for localized case than the delocalized 
case (figure 5(d)). Furthermore the presence of a symmetry weakens the repulsion at 
short distances but does not seem to affect the correlations at large distances (figures 
5(b),(c),(d)) (but as shown later, long range correlations are indeed symmetry-dependent). 
Our comparison of short-range correlations under various symmetry-conditions for many 
iV-values (61,99,149,199) seems to suggest the decrasing tendency of differences in the 
repulsion-strength under various parametric conditions with increasing N. We believe 
that no difference should survive in N — > 00 limit; this is also corroborated by our spacing- 
study given later. Moreover the short-range (r < 1) correlation between zeros is different 
from that for the two dimensional random distribution of points with R2(r) = 1 but 
behaves in a similar fashion for r > 1. Note here that the long range behaviour for all the 
cases approaches 1, although oscillating heavily around it, irrespective of the parameteric 
conditions or the type of the distribution. These oscillations are not insignificant; they 
should contain information about the differences in long-range correlations under various 
parametric conditions so as to agree with our number-variance study, given later. 

The i?2(r)-study suggests the short-range correlations to be free from the influence 
of various parameters. To understand this better, we carry out a study for the nearest 
neighbour spacing distribution P(s) of zeros (with all geometrical factors scaled out; see 
[10]) for both QKR as well as random functions. A comparison of the two reconfirms the 
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conclusion obtained by i? 2 (r)-study, that is, the short-range fluctuation-measures of zeros- 
distribution of random functions can very well model that of QKR (figures 6(a),(b),(c)). 
To understand the functional behaviour of -P(s), we fit the following equation 

P(s) = as a exp[-/3s b ] (14) 

where a and 6 are the fitting parameters and a and (3 are obtained by imposing following 

2 2 

normalization conditions, J Q P(s)ds = 1 and J Q sP(s)ds = 1. The fitting analysis carried 
out for four values of iV = 61, 99, 149, 199 leads us to expect a = 1.5 and 6 = 5. Using these 
values as a hypothesis for P(s), We perform the x 2 -analysis which gives us x 2 = 12.56 for 
N = 199 and 14 degrees of freedom at 5% level of significance, a value much less than the 
theoretical c/n 2 -value (= 23.68). This therefore indicates the correctness of our hypothesis. 
We also find that the calculated % 2 -value decreases with increasing iV (e.g. x 2 = 148.18 
for N = 99 and 14 degrees of freedom) which implies the tendency of P(s) to approach our 
hypothesis in N — > oo limit. 

Our P(s)-analysis also shows that the short range correlations between the zeros are 
not of the same type as those for a two-dimensional random distribution of points ( with 
P(s) as two dimensional Poisson distribution with a = and 6 = 2 [14]). In fact, it has 
been indicated [8] that the apperance of a term of the type of general n-body interaction 
in the dynamics of zeros should produce a strong short-distance correlations among zeros. 
Furthermore, for finite iV-cases, the repulsion between these zeros for small distances is 
stronger for localized case than the delocalized case. The presence of the repulsion leads 
zeros of the delocalized case to distribute uniformly in the phase space but the distribution 
is not a random one. In the localized case too, the existence of a stronger repulsion may 
be the cause of the pattern formation of zeros. Again the presence of a symmetry affects 
the spacing distribution for small-s values, by increasing the probability of small-spacings 
relative to no-symmetry case (figure 6(d)); this happens due to many very close-lying zeros 
on the symmetry- axis. But our analysis for various iV-values indicates the tendency of this 
probability to decrease with higher N, approaching the same behaviour as in no-symmetry 
case (i.e eq.(14)) in N — > oo limit. This is what we expect too as the fraction of zeros 
(iV( — 1/2)) on the symmetry-axis vanishes in the semiclassical regime. The fitting analysis 
followed by x 2 -analysis further confirms that, in N — > oo-limit, P(s) is not affected by the 
presence or absence of a symmetry, localization or derealization of quantum dynamics and 
can be described by eq.(14). (The calculated % 2 -value, if eq.(14) is taken as hypothesis, is 
12.56, 18.24 and 30.78 for 14 degrees of freedom of 1-symmetry, 2-symmetry and localized 
cases respectively. This is much less than the theoretical % 2 (= 23.68) at 5%level of 
significance for the first two cases, thus indicating correctness of our hypothesis. For the 
last case although the observed value is much greater than the theoretical one even for 
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N = 199 but again the decreasing % 2 with increasing N (e.g % 2 = 199.65 for N = 99 and 
14 degrees of freedom) suggests the validity of eq.(14) in N — > oo limit for this case too. 

The i?2(^)-study does not inform us much about the effect of various parameters on 
long range correlations. We study, therefore, another fluctuation measure, namely, the 
"number variance" n(r) which is defined as the variation in the number of zeros in an area 
of size r 2 , where r 2 is sufficiently large to include many zeros on an average and is two- 
dimensional analogue of the number variance quite often used for spectrum studies. For a 
random distribution of points distributed on a unit plane (with no correlation among them) , 
the variance increases as r 2 . If the variance remains smaller than r 2 as the mean number of 
zeros increases, the existence of long-range correlation in the spectrum is indicated. The 
figure 7. shows that the variance in the number of zeros for large r is increasing with 
r; the rate of increase, for both QKR and RF, is much slower than that of the random 
distribution of points. This reflects that the longe range correlations are indeed present 
among zeros of both QKR and RF, although stronger in the former as compared to latter, 
if the dynamics is localized or symmetries are present (figures 7(a),(b),(c)); the reverse is 
true if the dynamics is delocalized with no symmetry (figure 7(a)). As can be seen from 
figure 7(d), these correlations are strongest for delocalized dynamics with no symmetry, a 
little weaker for symmetry cases and weakest for localized case. A similar tendency was 
seen in case of the eigenvalues too [5,12]; the presence of a symmetry or localization of 
dynamics results in weakening of the long range correlations with a tendency to behave as 
in the case of Poission distribution. (This should not be surprising as Poission distribution 
results due to a large number of symmetries present in the system). 

In this paper, we have attempted to gain an insight in the complex world of wavefuc- 
tions by numerically analyzing the distribution of their zeros. Although the studies here 
have been carried out for one system only, we believe in general applicability of results ob- 
tained here. Our study indicates that the zeros for strongly chaotic cases tend to accquire 
uniform distribution in the phase space in the limit N — > oo; this distribution is different 
from that of the set of random points. The analysis of various fluctuation measures shows 
that, in the large A-limit, the short-range correlations of zeros of quantum chaotic systems 
are universal, parametric-independent and have the same nature as those for random func- 
tions case or random polynomials but long-range correlations seem to differ. We expect 
the nature of this deviation to be system-dependent; this expectation has its root in the 
earlier observation of similar tendency in the eigenvalues-case. Again, the localization vs 
derealization phenomena of quantum state reveals itself in the statistical behaviour of 
zeros too. The presence of a symmetry induces the changes in the distribution of zeros and 
also modifies the long range correlations. This again confirms that zeros contain a lot of 
information about the underlying classical as well as quantum dynamics and therefore can 
be used for a better understanding of the subject. But we have not as yet extracted all the 
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information about zeros. For example, we still have to understand as to how the motion 
of zeros is affected by the variation of various parameters of the system. We also need to 
analyze higher order correlation which can tell us more about their parametric and system 
dependence. It is also desirable to have a complete analytical formulation of statistical 
properties of these zeros as it will help us in a better understanding of wavefunctions. 
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Figure Captions 



Figure 1.: The superposition of the zeros of all eigenfunctions of U, 

(a) K = 20000, M = 1, N = 99, 7 = 0.7071, 9 = n/2N, 

(b) K = 10, M = 30, N = 99, 7 = 0.7071, # = n/2N. 

Figure 2.: The Husimi distribution of a single eigenfunction, 

(a) K = 10, M = 30, N = 61, 7 = 0.7071, # = n/2N, 

(b) K = 10, M = 30, JV = 61, 7 = 0.0, # = 0.0, 

(c) K = 10, M = 30, JV = 61, 7 = 0.0, 6 = 0.0, 

(d) K = 20000, M = 1, N = 61, 7 = 0.7071, # = vr/2iV. 

Figure 3.: The density p(q,p) of zeros in the phase plane with respect to q (averaged over 
all values of p and measured in units of < q >) and p (averaged over all values of q and 
measured in units of < p >), 

(a) & (b) K = 20000, M = 1,N = 99, 7 = 0.7071, O = 7r/2iV, (solid line), K = 
10, M = 30, iV = 99, 7 = 0.7071, 6 = tv/2N, (dashed line), RFC case (dotted line), 

(c) & (d) K = 20000, M = 1,N = 99, 7 = 0.0, O = 7r/2iV, (solid line), RFR case 
(dotted line), 



Figure 4.: The superposition of the zeros of all eigenfunctions of U, 

(a) K = 10, M = 30, iV = 61, 7 = 0.0, 9 = 0.0, 

(b) K = 10, M = 30, iV = 61, 7 = 0.0, 9 = n/2N. 

Figure 5. The histogram for Two-point correlation function R2 of the zeros (with unit 
mean spacing), 

(a) Delocalized without any symmetry (solid histogram) and RFC (dotted histogram), 

(b) Delocalized with one symmetry (solid histogram) and RFR (dotted histogram), 

(c) Delocalized with two symmetries (solid histogram) and RFA (dotted histogram), 

(d) with respect to small values of r, no symmetry case (solid line), 1-symmetry (dot 
line), 2-symmetry (small-dash line), localized case (big-dash line). Note here the difference 
between various curves is bigger than the finite-smaple error (ks 5 x 10~ 3 ) associated with 
each case. 



Figure 6. The Histogram for Nearest- Neighbour spacing distribution, (with unit mean 
spacing) 
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(a) Delocalized without any symmetry (solid histogram) and RFC (dotted histogram), 

(b) Delocalized with one symmetry (solid histogram) and RFR (dotted dotted his- 
togram) , 

(c) Delocalized with two symmetries (solid histogram) and RFA (dotted histogram), 

(d) Localized case with no symmetry (solid histogram) and Delocalized case with no 
symmetry (dotted histogram), 

Here solid curve in (a),(b),(c) is the fitting given by eq.(14). 

Figure 7. Number Variance n(r), 

(a) Delocalized without any symmetry (solid curve), Localized case (dashed curve) 
and RFC (dotted curve) , 

(b) Delocalized with one symmetry (solid curve) and RFR (dotted curve), 

(c) Delocalized with two symmetries (solid curve) and RFA (dotted curve). 
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